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We present results on the nucleon scalar, axial and tensor charges as well as on the momentum 
fraction, and the helicity and transversity moments. The pion momentum fraction is also presented. 

The computation of these key observables is carried out using lattice QCD simulations at a phys¬ 
ical value of the pion mass. The evaluation is based on gauge conhgurations generated with two 
degenerate sea quarks of twisted mass fermions with a clover term. We investigate excited states 
contributions with the nucleon quantum numbers by analyzing three sink-source time separations. 

We find that, for the scalar charge, excited states contribute significantly and to a less degree to the 
nucleon momentum fraction and helicity moment. Our result for the nucleon axial charge agrees 
with the experimental value. Furthermore, we predict a value of 1.027(62) in the MS scheme at 2 
GeV for the isovector nucleon tensor charge directly at the physical point. The pion momentum 
fraction is found to be {x)Z,-d ~ 0.214(15)(t9^) in the MS at 2 GeV. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw, 12.38.-t, 14.70.Dj 


I. INTRODUCTION 

The nucleon axial-vector coupling or nucleon axial 
charge gA is experimentally a well known quantity de¬ 
termined from the /^-decay of the neutron. It is a key 
parameter for understanding the chiral structure of the 
nucleon and a quantity that has been studied extensively 
in chiral effective theories mm- A description of baryon 
properties in chiral effective theory requires as an input 
gA and thus its value at the chiral limit and its depen¬ 
dence on the pion mass constitute important information 
that lattice QGD can provide. Its importance for phe¬ 
nomenology as well as the fact that it is rather straight¬ 
forward to compute in lattice QCD have made it one 
of the most studied quantities within different fermion 
discretization schemes mm- In lattice QCD, gA is deter¬ 
mined directly from the zero momentum transfer nucleon 
matrix element of the axial-vector current without requir¬ 
ing any extrapolation from finite momentum transfer cal¬ 


culations as, for example, is required for the anomalous 
magnetic moment of the nucleon. In addition, being an 
isovector quantity, it does not receive any contributions 
from the coupling of the current to closed quark loops 
and thus one only needs to compute the connected con¬ 
tribution with well established lattice QCD techniques. 
Therefore, gA has come to be regarded as a prime bench¬ 
mark quantity for the computation of lattice QCD matrix 
elements. Postdiction of the value of gA within lattice 
QCD is, therefore, regarded as an essential step before 
the reliable prediction of other couplings and form fac¬ 
tors for which the same formalism is used. 

Unlike gA, the nucleon scalar and tensor charges are 
not well known. Limits on the value of the scalar and 
tensor coupling constants arise from 0®" -T O®' nuclear 
decays and the radiative pion decay tt —)■ ez/y, respec¬ 
tively. They have become the focus of planned experi¬ 
ments to search for physics beyond the familiar weak in¬ 
teractions of the Standard Model sought in the decay of 
ultra-cold neutrons [TO]. The computation of the tensor 
charge is particularly timely since new experiments using 
polarized ®He/Proton at Jefferson lab aim at increasing 
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the experimental accuracy of its measurement by an or¬ 
der of magnitude m- In addition, experiments at LHC 
are expected to increase the limits to contributions aris¬ 
ing from tensor and scalar interactions by an order of 
magnitude making these observables interesting probes 
of new physics originating at the TeV scale. Computing 
the scalar charge will also provide input for dark matter 
searches. Experiments, which aim at a direct detection 
of dark matter, are based on measuring the recoil energy 
of a nucleon hit by a dark matter candidate. In many 
supersymmetric scenarios m and in some Kaluza-Klein 
extensions of the standard model mill] the dark matter 
nucleon interaction is mediated through a Higgs boson. 
In such a case the theoretical expression of the spin inde¬ 
pendent scattering amplitude at zero momentum transfer 
involves the quark content of the nucleon or the nucleon 
sigma-term, which is closely related to the scalar charge. 
In fact, this contributes the largest uncertainty on the 
nucleon dark matter cross section. Therefore, comput¬ 
ing the scalar gs and tensor gx charges of the nucleon 
within lattice QCD will provide useful input for the ongo¬ 
ing experimental searches for beyond the standard model 
physics. 


Another experimental frontier that provides informa¬ 
tion on the quark and gluon structure of a hadron, is the 
measurement of parton distribution functions (PDFs) in 
a variety of high energy processes such as deep-inelastic 
lepton scattering and Drell-Yan in hadron-hadron col¬ 
lisions. PDFs give, to leading twist, the probability 
of finding a specific parton in the hadron carrying cer¬ 
tain momentum and spin, in the infinite momentum 
frame. Their universal nature relies on factorization the¬ 
orems that allow differential cross-sections to be written 
in terms of a convolution of certain process-dependent 
coefficients that encode the hard perturbative physics 
and process-independent PDFs that describe the soft, 
non-perturbative physics at a factorization energy scale 
g, [13[TB]. Because these PDFs are light-cone correlation 
functions it is not straightforward to calculate them di¬ 
rectly in Euclidean space. Instead, one calculates Mellin 
moments of the PDFs expressed in terms of hadron ma¬ 
trix elements of local operators, which through the oper¬ 
ator product expansion are related to the original light- 
cone correlation functions. Mellin moments are measured 
or extracted from phenomenological analyses in deep- 
inelastic scattering experiments and thus they can be 
directly compared to lattice results when converted to 
the same energy scale g. 


{X)q 

{x)Aq 

{x)Sq 



[q{x) + q{x)] dx 



[Aq{x) — Aq{x)] dx 



[Sq{x) + Sq{x)] dx , 


( 1 ) 

( 2 ) 

(3) 


where qi and q^ correspond respectively, to quarks with 
helicity aligned and anti-aligned with that of a longitu¬ 
dinally polarized target, and q-f and q± correspond to 
quarks with spin aligned and anti-aligned with that of a 
transversely polarized target. These moments, at lead¬ 
ing twist, can be extracted from the hadron matrix el¬ 
ements of one-derivative vector, axial-vector and tensor 
operators at zero momentum transfer. Thus, they consti¬ 
tute the next level of observables in terms of complexity 
that can be computed in lattice QCD after the coupling 
constants that do not involve derivative operators. The 
unpolarized and polarized moments {x)q and {x)Aq of 
the nucleon are measured experimentally and thus lat¬ 
tice QCD provides a postdiction, while a computation of 
the nucleon transversity {x)sq provides a prediction. It is 
worth mentioning a new approach proposed recently for 
measuring directly the PDFs within lattice QCD nanH], 
which is currently under investigation mm- 

In this paper, we extend the analysis of meson masses, 
the muon anomalous magnetic moment g — 2 and the me¬ 
son decay constants considered in Ref. |3T] , to the nucleon 
matrix elements for the three first Mellin moments while 
for the pion we compute the momentum fraction. While 
the present paper builds on the methodology developed 
in Refs. this work presents the first evaluation 

of these six quantities directly at the physical value of 
the pion mass. This is a substantial step forward since it 
avoids chiral extrapolations, which are often difficult and 
can lead to rather large systematic uncertainties. 

The paper is organized as follows: In section H we 
define the nucleon and pion matrix elements, in section 
HI we explain the lattice methodology, in section IV we 
give the simulation details and in section V our results. 
Section VI summarizes our findings and gives our con¬ 
clusions. 


II. MATRIX ELEMENTS 

A. Nucleon 


In this work, we consider the three first moments that 
one can construct, namely the first moment of the spin- 
independent (or unpolarized) q = qi + q^, helicity (or 
polarized) Aq = qi — q^, and transversity Sq = q-f + q± 
distributions, which are define as follows: 


We are interested in extracting the forward nucleon 
matrix elements {N{p)\0\N(p)), with p the nucleon ini¬ 
tial and final momentum. We consider the complete set 
of local and one-derivative operators, yielding a non-zero 
result. The local scalar, axial-vector, and tensor opera¬ 
tors are: 
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Os’^=q^q, C>Aa = (?757'"y9 , = qa^'^^q . 

(4) 

We do not consider the vector operator 'tlj{x)jmUip{x) 
since this yields the renormalization constant Zy, which 
we calculate separately using our RI-MOM setup, as ex¬ 
plained in section V. If one instead uses the lattice con¬ 
served Noether current, which we typically do in in our 
computation of the nucleon electromagnetic form factors, 
we the forward matrix element trivially yields the electric 
charge. 

The one-derivative vector, axial-vector, and tensor op¬ 
erators are given by 

OfiT. = 

0^7 = (5) 

where q = {u, d), 

0/1= 2 i^Dfi — = 2 4" V/j) 

and Vm (Vp is the usual forward (backward) derivative 
on the lattice. The curly (square) brackets represent a 
symmetrization (anti-symmetrization) over pairs of in¬ 
dices, with the symmetrization accompanied by subtrac¬ 
tion of the trace. 


In what follows all expressions will be given in Eu¬ 
clidean time. For example, the one-derivative vector cur¬ 
rent that will be used in the computations of both nu¬ 
cleon and pion momentum fractions, in Euclidean time 
and setting fj, = v = A, is given by: 



In this work, we consider the isovector quantities ob¬ 
tained from Eqs. Q and ([^ by using the Pauli matrix 
r^. We also consider the isoscalar combination obtained 
by replacing r“ by unity. The individual up- and down- 
quark combinations can be extracted form the isovector 
and isoscalar quantities which are equivalent to replacing 
t“ with the projectors onto the up- or down- quarks. The 
isoscalar combination and the up- and down-quark con¬ 
tributions receive disconnected contributions. Our high 
statistics study using a.n Nf = 2-I-1-I-1 ensemble of 
twisted mass fermions with pion mass of 373 MeV has 
shown that the disconnected contributions for the tensor 
isoscalar charge and the isoscalar first moments are very 
small compared to the connected HZlllHl- In the same 
study, the disconnected contributions to the isoscalar ax¬ 
ial and scalar charge were found to be about (7- I0)% of 
the connected. In this work, we will only compute the 
connected contributions. The disconnected contributions 
will need at least an order of magnitude more statistics 
and will be presented in a follow-up publication. 

For zero momentum transfer, the nucleon matrix ele¬ 
ments of the local operators in Eq. Q can be decomposed 
in the following form factors: 


(iV(p,s')|Os|iV(p,s)) 

{N{p,s')mN{p,s)) 

{N{p,s')\Oli!'\N{p,s)) 


un { p , s ') -Gs(0) un { p , s ), 
iuN{p,s') un{p,s), 

un { p , s ') ^ATio(0)cr'^'' un { p , s ). 


(7) 

( 8 ) 
(9) 


Thus, the scalar matrix element at zero momentum 
transfer yields the form factor GsiO) = gs^ the local 
axial-vector Gyi(O) = qa and the local tensor matrix el¬ 
ement yields Ano(O) = gr- In all these quantities the 
operators are either the isovector or isoscalar combina¬ 
tions, or individual up- or down-quark contributions. At 
non-zero momentum, additional form-factors arise in the 
decomposition of Eqs. Q and Q. Namely, the induced 
pseudo-scalar Gp{Q‘^) appears as the second form factor 


in the decomposition of the matrix element of the axial- 
vector and the form factors Stio(Q^) and Axio{Q^) ap¬ 
pear in the decomposition of the nucleon matrix element 
of the tensor operator, where is the momentum trans¬ 
fer square in Euclidean time. These cannot be extracted 
at zero momentum transfer and will not be considered in 
this work. 

The corresponding decomposition for the one- 
derivative operators in Eq. ([^ is given by: 
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{Nip,s')\0^’'\Nip,s)) 

{N{p,s')\0^/\N{p,s)) 

{N{p,s')\0!i!'P\N{p,s)) 


un{p,s') -A2 o( 0)7 ^'"p^^ un{p,s), 
iuN{p,s') ^^20(0)7^'"^'"^® un{p,s), 
iUN{p,s') ^AT 2 o{ 0 )cr^''^''^P^^ un{p,s). 


( 10 ) 

( 11 ) 

( 12 ) 


The momentum fraction, helicity moment, and the 
transversity moment are obtained from the above for¬ 
ward matrix elements by {x)q = ^ 2 o(l*)> {^)^q — ^ 2 o( 0 )j 
and {x)sq = ^^ 20 ( 0 ) respectively. Here we use the 
generic symbol q to denote the quark combination, where 
q = u + d will denote the isoscalar combination, q = u — d 
will denote the isovector combination and q = it or 
q = d denotes the individual up- and down-quark con¬ 
tributions. For instance, the isovector helicity moment 
will be denoted as {x)Au-Ad = A2 q'^{0). For uniformity 
in our notation we will also write for the nucleon 

axial charge, despite the fact that the measured axial 
charge is understood to be an isovector quantity. 

B. Pion 

The isovector momentum fraction of the pion {x)'!^^^, 
can be extracted from the corresponding pion matrix el¬ 
ement of the one-derivative vector operator. Specifically 
we use the following operator, sometimes also denoted as 


^v2bt 

0^\x) = lo^\{x) (13) 

where 0^4 (x) is given in Eq. As in the case of the 
nucleon, no external momentum is needed in our calcula¬ 
tion, which is advantageous since an external momentum 
increases the noise to signal ratio. 


III. LATTICE METHODOLOGY 

A. Correlation functions 

In order to compute hadron matrix elements we need 
to calculate the appropriate three-point function. We 
first present the setup for the nucleon matrix elements 
for the special case q = 0. The three-point function is 
then given by 


^ ^ e-hx.-xo). 


where Xq, Xins and Xg are the source, insertion and sink 
coordinates respectively. In order to cancel unknown 
overlaps of the interpolating field with the nucleon state 
as well as the time evolution in Euclidean time we con¬ 
struct ratios of the three-point function with the two- 
point function, which is given by 

G2pt(0,ts) = ^r^„(Ja(xs,ts)J/3(xo,to))- (15) 

Xs 

The projection matrices are 

r4 = i(i + 74), r'^ = r4z757,. (le) 

We use the proton interpolating operators: 

Ja{x) = e“'"=<(a;)[u^'’(x)C'75d‘=(x)] (17) 

with a, b and c denoting color components. We employ 
Gaussian smeared quark fields [30] to increase the 


^s)Or"’"'’'"" (Xins, tins)^/3(Xo, <o)) (14) 


overlap with the proton state and decrease overlap with 
excited states. The smeared interpolating fields are given 
by 

9“mear(^>x) = ^ (x,y; [/(t)) (t, y), (18) 

y 

F = + 

3 

H{x,y,u{t)) = '^[Ui{x)S,;^y_i + Uj{x-i)Sg;^y+i]. 

i -1 

We apply APE-smearing to the gauge fields Uy entering 
the hopping matrix H. The parameters for the Gaussian 
smearing ac and Nq are optimized using the nucleon 
ground state |3T| such as to give a root mean square 
radius of about 0.5 fm. We use {NQ,ac) = (50,4) and 
(YapejOape) = (50,0.5). 
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FIG. 1: Connected (upper) and disconnected (lower) contri¬ 
butions to nucleon three-point functions. 


For the case of isovector quantities, the so-called dis¬ 
connected contributions arising from the coupling of the 
operators to a sea quark shown schematically in the lower 
panel of Fig. are zero in the isospin limit up to lattice 
cut-off effects. Since we work with an automatic 0{a)- 
improved action we expect cut-off effects to be small for 
our action and lattice spacing. Thus, these correlators 
can be calculated by evaluating the connected diagram 
of Fig. (upper panel) for which we employ sequen¬ 
tial inversions through the sink [32]. For the case of 
isoscalar quantities, the disconnected diagrams do not 
vanish and need to be computed. The calculation of 
disconnected contributions needs special techniques and 
at least an order of magnitude more statistics than the 
connected ones |5Sj. In a high-statistics analysis using 
Nf = 2-I-I-I-1 twisted mass fermions at pion mass of 
373 MeV we found that they contribute about 7% to the 
isoscalar axial charge while they are negligible for the ten¬ 
sor charge, {x)q and {x)Aq [IZ]- In this work we restrict 
ourselves to the calculation of the connected contribu¬ 
tions. 

For evaluation of the connected contribution we take 
the nucleon creation operator at a position xq = ( xq , to), 
referred to as the source position, which is randomized 
over gauge-configurations to reduce autocorrelations, and 
the annihilation operator at a later time tg, the sink po¬ 
sition Xs- The current couples to a quark at an interme¬ 
diate time tins) the insertion time. In this calculation, we 
fix the sink and source to momentum p=0 and thus the 
current carries zero momentum. We employ the sequen¬ 


tial inversion method through the sink requiring one set 
of sequential inversions per choice of the sink time-slice 
tg and sink projector. Thus, within this approach, at 
a fixed sink-source time separation we obtain results for 
all insertion times as well as for any operator 
We perform separate inversions for the four projection 
matrices F^ and F^ given in Eq. (16). 


Using the two- and three-point functions of Eqs. (14) 
and (15) we form the ratio 


i7(F^,tg,ti„s) = 


G'3pt(r^, 0, tg, tins) 

G2pt(0) tg) 


(19) 


where all time separations are given relative to the source 
time to- Eor large time separations, the two-point func¬ 
tion in the denominator cancels unknown overlaps of the 
nucleon interpolating operators with the nucleon spinors 
as well as the Euclidean time evolution, such that the 
desired matrix element is isolated. However, in order to 
identify when the large time limit sets in, one has to care¬ 
fully study the time dependence of the ratio since excited 
states contamination can affect the value of the ratio. For 
arbitrary times the contributions from excited states in 
the ratio are 


-R(r^,tg,tins) oc 




( 20 ) 


where \n) is the eigenstate of the QCD Hamiltonian 
with the quantum numbers of the nucleon and is the 
energy of the state in the rest frame of the nucleon. De¬ 
noting with |0) = \N) the nucleon ground state, with 
|1) = \N') and |2) = \N") the first and second excited 
states, and with A = — rriN and A' = En" — Ejqi 


their respective energy gaps, the ratio in Eq. (20) yields 


) tg , tins ) CC 


I /n/'„-A'(ti, 


1 -b |C|2e-^(*.-‘o) -p o(e-A'(t.-io)) 


( 21 ) 


where A4 = (iV|(!lr|A) is the desired matrix element, 
C = and TZ = C{N'\Or\N). If the exponential 

terms are small compared to Ai and to unity, then we 
have what is referred to as ground-state dominance and 
the ratio yields the desired ground state matrix element. 

The computation of the pion momentum fraction is 
carried out along the same lines as for the nucleon. To 
extract the desired matrix element (7r(0)|(!l"‘'*|7r(0)) we 
construct the ratio of the appropriate three-point func¬ 
tion with the pion two-point function with a source at to 
and a sink at tg: 


R'^{tg,tins) = 


GitMts 


^ins) 


^2pt(^s) 


( 22 ) 
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with 


B. Ensuring ground state dominance 


G3pt(0:^s,tins) = jUito)) ■ 

(23) 

with J^+(x) = d[x)y^u{x) (J^-{x) = u{x)y^d{x)) is the 
interpolating field of 7r+ (tt”). As in the case of the nu¬ 
cleon, in the isovector combination up to lattice artifacts, 
the disconnected contributions vanish and will thus be 
dropped. 

The pion two- and three-point correlators, unlike those 
of the nucleon, are evaluated by using a stochastic 
time-slice source (Z(2)-noise in both real and imaginary 
part) [33U35] for all color, spin and spatial indices. This 
method, which is particularly suited for the pion, has 
been first applied to moments of parton distribution func¬ 
tions in Ref. [36] . The quark propagator 5'^(y) is ob¬ 
tained by solving 


To extract the nucleon matrix element from the ra¬ 


tio defined in Eq. (19) one needs to make sure that the 


contribution of the terms due to the excited states in nu¬ 


merator and denominator of Eq. (21), the so called con¬ 


tamination due to excited states, is negligible. We will 
employ three methods to check for ground state domi¬ 
nance, as described below. 

In the first method, which we will refer to as the plateau 
method, one probes the region for which A(ts — tins) ^ 1 
and A (tins —to) ^ 1 such that excited state contributions 
are much smaller than the contribution of the ground 
state. Within this time interval the ratio becomes time 
independent and the time range where this happens is 
referred to as the plateau region. Fitting the ratio 

R{T\ts,tins) n(r^) (27) 


y) Sl{y) = C(z)“ (source at <o) (24) 

y 

for 5. ^(z)^ is a Z{2) random source satisfying 

(^(x)“ e(y)^). = o, 

where (•)n denotes the average over many random 
sources. Using S we can define a so-called sequential 
or generalized propagator ^'a{y) from |37| 


XI = l5S'^{z)5^n,ts (sink at ts). 

y 

(26) 

This method represents a generalization of the one-end- 
trick [38] to moments of parton distribution functions. 
Its clear advantage is an increased signal to noise ratio at 
reduced computational costs at least when it is applied 
for meson observables. With a point source, 24 inver¬ 
sions per gauge configuration are needed: 12 (3 colors x 
4 spins) for the quark propagator and 12 for the gener¬ 
alized propagator. With the stochastic source discussed 
above, only two inversions are needed: one for the quark 
propagator and another one for the generalized propaga¬ 
tor. For a comparison of stochastic versus point sources 
we refer to Ref. [36] . 

The stochastic method described above can be adapted 
to work for other mesons. However, for moments of nu¬ 
cleon parton distribution functions we found no improve¬ 
ment in the signal to noise ratio for a comparable com¬ 
putational effort. 

For {x)^^ it is sufficient to fix tg — to = T/2, where 
T is the temporal extent of the lattice. As in the case of 
the nucleon, the value of to is chosen randomly on every 
gauge configuration in order to reduce autocorrelation. 


over tins within this plateau region one obtains the 
plateau value, which is the desired matrix element M.. 
To ensure excited state suppression one repeats this pro¬ 
cedure for multiple values of tg , checking that the plateau 
value does not change. However, the statistical errors 
grow exponentially with tg , which means that as the sink- 
source time separation increases the signal is lost as com¬ 
pared to the statistical noise making it difficult to detect 
any time-dependence. Increasing tg therefore requires a 
corresponding increase in statistics if this check is to be 
useful [25]. 

The second approach is to use the summation method 
proposed some time ago ^S] and recently applied to the 
study of the nucleon axial charge m- One sums the 
ratio over the time of the insertion. 


ts—T 

^sum(rA)^ ^ R{r\tg,tins), (28) 

tinB=*0-|-T 

with T selected such that contact terms are not included, 
i.e. T = 1 for local operators and r = 2 for derivative 
operators. The sum over the excited state contributions 
given in Eq. © is a geometric series and can easily be 
summed to yield 

^sum(rA) (29) 


with C a constant independent of tg. The advantage over 
the plateau method is that excited state contamination 
is suppressed by a larger factor [A(ts — to) as opposed to 
A{tg — tins) or A(tins —io)]- However, the extraction of M 
requires a fit to two parameters, resulting in general in 
larger statistical uncertainties. Nevertheless, this method 
provides a good consistency check of our results. 

A third approach to extract the desired matrix el¬ 
ement is to take into account in the fit the contribu¬ 
tion of the first excited state in Eq. (21). In this case, 
we simultaneously fit the two- and three-point correla¬ 
tion functions obtained from the lattice including the 
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ground state and the first excited state contributions. 
This is done by performing a combined fit to all sink- 
source separations and to both correlation functions with 
tins and ts as independent variables. Like for the sum¬ 
mation method, we exclude the contact terms, i.e. for 
tins S [to-l-1, ts —1] for the scalar, axial and tensor charges, 
and tins € [tg -|- 2, ts — 2] for the momentum fraction, po¬ 
larized moment and trasversity moment, which include a 
derivative. We will refer to this method as the two-state 
fit method. 

In this work, we consider agreement among the above 
three methods yielding the same value for M. as our cri¬ 
terion that excited states are sufficiently damped out. 

If one has ground state dominance the nucleon matrix 
elements of the scalar, axial and tensor local operators, at 
zero momentum transfer and Euclidean time are related 
to the ratio as follows: 

n5(r") = f 

n^^(rfc) = 

n^(rfc) = (30) 

The corresponding expressions for the vector, axial and 
tensor one-derivative operators are: 

nt4(r4) = 

= 

ni"(rfc) = 

= 

(31) 

Note that after symmetrization and subtraction of the 
trace as indicated in Eq. ([^, only one of the two expres¬ 
sions for {x)u±d is independent. 

Eor the case of the pion we only present We 

consider the largest sink-source separation possible on 
each lattice, namely ts — to = T /2. This is possible for 
the pion since its two point function has constant signal 
to noise ratio independently of ts — to. Therefore, we 
extract the pion momentum fraction using the plateau 
method at this single value of the sink-source separation: 

R^{ts,Ur,s) n" ( 32 ) 

A(tins — 

where we use A to generically denote the energy gap be¬ 
tween the energy of the first excited state and the ground 
state of the hadron of interest in its rest frame. Given 
ground state dominance, the pion momentum fraction is 
at zero momentum transfer and Eucledian time obtained 
from the ratio via: 

= ( 33 ) 


^ {x)u±d 


niN 

^ Wu±d 


2 (x) Au± Ad 


IV. SIMULATION DETAILS 

We use the (maximally) twisted mass fermion (TME) 
formulation of lattice QCD [41], which is particularly 
suited for hadron structure calculations since it provides 
automatic 0{a) improvement requiring no operator mod¬ 
ification [421445] . Twisted mass ensembles with two de¬ 
generate flavors of light sea quarks {Nf = 2) as well as 
ensembles including the strange and charm sea quarks 
{Nf = 2 4-1-1-1) are produced by the European Twisted 
Mass Collaboration (ETMC) and technical details on the 
simulations can be found in Refs. [15IHH] and [JS] respec¬ 
tively. This work focuses on the analysis of gauge config¬ 
urations produced using two degenerate flavors of twisted 
mass light sea quarks {Nf = 2) including a clover term. 
Eor the gauge action we use the Iwasaki action. The pa¬ 
rameters of the four ensembles considered in this work 
are given in Table |l| More details on the choice of action 
and the simulations are given in Refs. [snuoiiii]. 


TABLE I: Input parameters of our new lattice ensembles used 
in this work. For each ensemble we give the lattice size, the 
bare quark mass (a/i) and the corresponding pion mass (m.^-). 
These ensembles use TMF at one value of f with a clover 
term with csw ~ 1.57551. The lattice spacing given in the 
table is determined using the nucleon mass as explained in 
the text. 


P = 2.1, a = 0.093(1) fm, rp/a = 5.32(5) 

0.006 0.003 

V) 0.338(9) 0.244( 

0.006 

V) 0.335(9) 

0.0009 

(GeV) 0.1312(13) 


24^ 

X 

00 

= 2.23 

fm 

a/i 






m^r 

32^ 

X 

64, L 

= 2.97 

fm 

a/i 






m^r 

48^ 

X 

96, L - 

II 

fm 

a/i 






rUn 


For the nucleon structure observables, we analyze the 
ensemble with afj, = 0.0009. We will refer to this ensem¬ 
ble as the physical ensemble and speak in what follows of 
the physical point. For the case of the pion momentum 
fraction, we use all four ensembles of TMF with a clover 
term. 

Although the observables of interest in this work are 
dimensionless and do not depend on the lattice spacing, 
it is useful to study their dependence on the pion mass, 
which is a dimensionfull quantity. In Ref. |2I| . the lat¬ 
tice spacing for the new Nf = 2 ensembles with the clover 
term was determined using gluonic quantities as well as 
the pion and kaon decay constants. Another determina¬ 
tion of the lattice spacing mentioned in Ref. [2T] is via 
the nucleon mass. 

The physical value of the nucleon mass was used as 
an input for the determination of the lattice spacings 
in our previous analysis of the Nf = 2 [22j [^ and 
Nf = 2-\-l-\-l [54] twisted mass ensembles. Each set of 
Nf = 2-I-I-I-1 and Nf = 2 ensembles involved three values 
of the lattice spacing. Since those simulations involved 










TABLE II: The lattice spacing and value of the scale parame¬ 
ter, ro |52], for the Nf = 2 and Nf = 2-|-1-1-1 TMF ensembles 
as well as for the new Nf = 2 TMF ensembles with the clover 
term, determined using the nucleon mass as explained in the 
text. The first error is statistical. The second error is the 
difference in the value when discarding ensembles with pion 
mass larger than 300 MeV. 

Nf = 2 


P 

a (fm) 
ro (fm) 


3.9 4.05 4.2 

0.088(2)(2) 0.071(2)(1) 0.056(2)(1) 
0.458(10)(1) 0.467(12)(7) 0.465(13)(7) 


Nf = 2 + 1 + 1 


P 

a (fm) 
ro (fm) 


1.90 

0.094(1)(2) 

0.501(7)(9) 


1.95 

0.082(1)(2) 

0.492(6)(3) 


2.1 

0.065(1)(1) 

0.499(6)(5) 


Nf = 2 with csw = 1.57551 


P 

a (fm) 
ro (fm) 


2.1 

0.093(1)(0) 

0.493(5)(0) 



[GeV^] 


FIG. 2: The ratio of the nucleon mass to the pion mass as a 
function of the pion mass squared. For determining the pion 
mass squared the scale is set using the nucleon mass at the 
physical point as described in the text. The fit only used the 
Nf = 2 + 1 + 1 ensembles without a clover term (filled circles, 
diamonds and squares). The plot also shows the Nf = 2 
TMF results (open circles, diamonds and squares) and the 
Nf = 2 ensemble with a clover term at the physical point 
(filled triangle). For the latter action {Nf = 2 with a clover 
term) we restrict our analysis of nucleon observables to the 
ensemble simulated at a physical value of the pion mass only. 


larger than physical light quark masses a chiral extrapola¬ 
tion was needed. We used the lowest order heavy baryon 
chiral perturbation theory expression, given by [55] 


arriN = am% — 4(ci/a)(aTO.n.)^ — 


35i 


167r(a/^)2 


{am^f 


(34) 

which is well-established within baryon chiral perturba¬ 


tion theory, is the value of the nucleon mass in the 
chiral limit and —4ci gives the cr-term written in units of 
the lattice spacing. The fit was constrained to reproduce 
the physical nucleon mass, by fixing the value of ci. In¬ 
cluding an a^-term in Eq. (34) had a negligible effect on 


the fit showing that indeed cut-off effects are small [54] 
for lattice spacings smaller than 0.1 fm. This justified 
the utilization of continuum chiral perturbation theory 
to determine the three lattice spacings for each Nf = 2 
oriV/ = 2 + l + lby simultaneously htting each set of 
17 Nf = 2 + 1 + 1 or 11 Nf = 2 ensembles. The val¬ 
ues of the nucleon mass used are taken from Ref. [^ for 
the Nf =2 + 1 + 1 ensembles and from Ref. (32) for the 
Nf = 2 ensembles. These values of the lattice spacings 
are used to obtain the pion mass for these ensembles. 

For the physical ensemble using the values 


= 0.06196(9) arriN = 0.440(4), 


(35) 


and assuming that we are exactly at the physical point 
we find a = 0.0925(8) fm where the average nucleon mass 
rriN = 0.939 GeV is used as an input. With this lattice 
spacing we find TOtt = 0.1323(12) MeV, where the largest 
part of the error comes from the error on the lattice spac¬ 
ing. This is about 5% less than the average physical pion 
mass. Using the values of Eq. (35) we find for the ra¬ 


tio mjv/m,r± = 7.10(6) compared to the physical value 
of 0.939/0.138 = 6.8, which again differs by less than 5% 
from the physical value. In order to check what the effect 
of a possible small mismatch in the pion mass would be 
on the lattice spacing, we use the fit extracted from the 
Nf = 2 + 1 + 1 ensembles to interpolate to the physical 
value of pion mass. This is done by making a combined 
fit of the 17 Nf = 2 + 1 + 1 ensembles with their three 
lattice spacings and the lattice spacing for the physical 
ensemble as well as as fit parameters. The fit yields 
X^/d.o.f = 1.6 for d.o.f = 12, which is a reasonable value. 
We find a value of a = 0.093(1) fm for the physical en¬ 
semble, consistent with the determination using Eq. (351, 


while the lattice spacings for the Nf = 2 + 1 + 1 remain 
unchanged compared to the values obtained when the 
physical ensemble was not included. Using a = 0.093(1) 
we find to^± = 0.1312(13) GeV, which is consistent with 
the value extracted from Eq. Excluding from the 
fit pion masses larger than 300 MeV yields consistent 
results for the lattice spacings of the Nf = 2 + 1 + 1 
ensembles while it does not change the value of the lat¬ 
tice spacing at the physical point. We note that if we 
fit using the Nf = 2 ensembles |22l |53] instead of the 
Nf = 2 + 1 + 1 ensembles the value of a = 0.093(1) fm is 
unchanged. This indicates that the mild interpolation is 
very robust. In Fig. we show the ratio of the nucleon to 
pion mass TOAr/m,n.±, which is a dimensionless observable 
determined purely from lattice quantities. We note that 
the values of the lattice spacings affect only the determi¬ 
nation of the pion mass plotted as the x-axis. The curve 
shown in Fig.j^is the fit to the ratio performed on the 17 
^^■ = 2 + 1 + 1 ensembles alone. The resulting chiral fit 
using m,r < 500 MeV yields y^/d.o.f = 1.4 and describes 
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very well the data. In the figure we also include the val¬ 
ues for the ratio for the Nf = 2 ensembles, which also 
fall on the same curve. This can be taken as an indica¬ 
tion that indeed strange and charm sea quark effects are 
small for the nucleon sector. The consistency of our new 
result is demonstrated in Fig. by the fact that the ra¬ 
tio mjv/m 7 r± for our physical ensemble falls on the curve 
determined from fitting the Nf = 2-|-l-l-l alone. Fig. 
provides a nice demonstration of the negligible effect M 
lattice artifacts on the m 7 v/m.,r± ratio. 

We note that the value of the lattice spacing deter¬ 
mined from the nucleon mass analysis is fully consistent 
with the one determined from gluonic quantities such 
as the one related to the static quark-antiquark poten¬ 
tial, To, and the ones related to the action density renor¬ 
malised through the gradient flow. It is, however, larger 
by about 1% as compared to that extracted using [21]. 
This was also observed in our analysis of Nf = 2 and 
Nf = 2-I-1-I-1 TMF ensembles [S5]. In Table we 
collect the lattice spacings for all the TMF ensembles de¬ 
termined using the nucleon mass and Eq. (34). We take 
as a systematic error due to the chiral extrapolation the 
shift in the mean value when discarding ensembles with 
pion mass greater than 300 MeV. For completeness we 
also give the values of rg determined from the nucleon 
mass in the same way as the lattice spacings, although 
they are not needed in this work. In what follows we 
use, for the physical ensemble, the value a = 0.093(1). 
The lattice spacings given in Table [Tl| are used to convert 
the pion mass to physical units. No other physical quan¬ 
tity presented in this work is affected by the value of the 
lattice spacings. 


V. RESULTS 

For the nucleon observables we analyze 96 gauge field 
configurations with 16 randomly chosen positions for 
each configuration yielding a total of 1536 measurements. 
For the nucleon observables, we use three sink-source 
separations for both the plateau and the summation 
methods, namely t^/a =10, 12, and 14, correspond¬ 
ing to approximately 0.9 fm, 1.1 fm, and 1.3 fm. For 
all separations we have 1536 measurements, by comput¬ 
ing the required two- and three-point correlation func¬ 
tions. First results on these quantities were presented in 
Refs. [571158] • For the pion we use the largest possible 
time separation namely T/2. 


A. Renormalization 

We determine the renormalization functions for the 
lattice matrix elements non-perturbatively, in the RF- 
MOM scheme employing a momentum source [59]. For 
the computation of the renormalization functions of the 
Nf = 2 + 1 + 1 ensembles we employed Nf = 4 simula¬ 


tions for at least three different values of the pion mass 
taking the chiral limit. A similar analysis was performed 
for the Nf = 2 TMF ensembles as well as for our new 
Nf = 2 TMF ensembles that include the clover term us¬ 
ing the ensembles with afi = 0.006,0.003 and 0.0009, the 
latter being at the physical pion point. In Refs. (60] |6T] 
we carried out a perturbative subtraction of 0{a^) terms 
that subtracts the leading cut-off effects yielding only a 
very weak dependence of the renormalization factors on 
{apY for which the (ap)^ —>■ 0 limit can be reliably taken. 
In this work, we reduce even further the 0{a?) contribu¬ 
tions by subtracting lattice artifacts computed perturba- 
tively to one-loop and to all orders in the lattice spacing, 
0{g'^ a°°), so that we eliminate a large part of the cut-off 
effects. In fig. [^ we show the results on the axial and ten¬ 
sor renormalization functions after subtraction. As can 
be seen, lattice artifacts are practially removed allowing a 
robust extrapolation to (ap)^ = 0. Due to the good qual¬ 
ity of the plateaus after the subtraction of the 0{g'^ a°°) 
any choice for the fit within the non-perturbative region 
(ap)^ £ (2 — 7) yields consistent results. Details on this 
computation can be found in Ref. [62] . 



FIG. 3: Results for the axial (upper) and tensor(lower) renor¬ 
malization functions as a function of the momentum square 
in lattice units. The open (black) circles are the unsubtracted 
results while the filled triangles (magenta) show the data af¬ 
ter 0{g^ a°°)-terms are subtracted. The filled (red) circle at 
(ap)^ = 0 is the value extracted by fitting to the plateau 
region [2-7] the subtracted data. 


Our previous chiral extrapolations have shown that for 
all renormalization functions except Zp the pion mass 
dependence is very weak. For the physical ensemble 
we compute Zp for three pion masses corresponding to 
ap = 0.006,0.003 and 0.0009 and performed the pole 
subtraction. Our value is given in Table III [S2]. The 


scheme and scale dependent renormalization functions 
are converted in the MS-scheme at a scale of 2 GeV, 
using the intermediate RGI scheme. 

We collect the values of all relevant renormalization 
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functions in Table III converting the scale dependent 
renormalization function in the MS scheme at a scale 
= 2 GeV, which is applicable to all except Z^. The 
systematic error is computed by varying the interval for 
the continuum extrapolation (ap)^ —>■ 0. The values of 
Zp for the fVy = 2 + 1 + 1 ensembles are taken from 
Refs. [561163] where the pole subtraction was performed, 


while for the new ensembles with the clover term we use 


the result of this paper. The values given in Table III are 
used to renormalize the lattice matrix elements studied 
in this work. More details are reported in Ref. [52] ■ For 
the Nf = 2 TMF ensembles without the clover term we 
do not calculate the scalar charge and transversity and 
therefore the renormalization functions are not given. 


TABLE III: Renormalization fnnctions for the ensembles used in this work. They are given in the twisted basis. They are the 
same in the physical basis except for Zp, which renormalizes the scalar operator in the physical basis. The renormalization 
functions for the local axial-vector, scalar and tensor operators are given in columns two, three and four, respectively. The 
three last columns give the renormalization functions for the derivative vector, axial vector and tensoroperators. The first error 
is statistical and the second error systematic. 



Za 

/7MS 

Zp 

/7MS ^MS ^MS 

Zt Zp,v ^DA 

^MS 

Zdt 

CSI 

II 

3.90 

4.05 

4.20 

0.769(2)(1) 

0.787(1)(1) 

0.791(1)(1) 


0.758(2)(4) 1.028(2)(6) 1.102(5)(7) 
0.796(1)(3) 1.080(2)(11) 1.161(4)(13) 
0.814(1)(3) 1.087(3)(12) 1.164(3)(6) 





II 


1.90 0.7474(6) (4) 
1.95 0.7556(5)(85) 
2.10 0.7744(7)(31) 

0.529(7)(45) 

0.509(4)(37) 

0.516(2)(29) 

0.7154(6)(6) 1.0268(26)(103) 1.1170(54)(223) 1.0965(90)(278) 
0.7483(6)(94) 1.0624(108)(33) 1.1555(36)(289) 1.1727(121)(73) 
0.7875(9)(15) 1.0991(29)(55) 1.1819(47)(147) 1.1822(59)(118) 

Nf=2 + csw 

2.10 

0.7910(4)(5) 

0.5012(75)(258) 0.8551(2)(15) 1.1251(27)(17) 1.1357(20)(205) 1.1472(121)(48) 


The renormalization functions are given for the twisted 
basis. Going from the twisted to the physical basis affects 
only the renormalization function for the scalar charge, 
which, in the twisted basis, is renormalized with Zp. Fur¬ 
thermore, since disconnected contributions are neglected, 
the isovector and isoscalar are renormalized using the 
same renormalization functions. All our results on the 
scalar and tensor charges and on the moments of PDFs 
are given in the MS scheme at an energy scale of 2 GeV. 


B. Nucleon scalar, axial and tensor charges 

In what follows we will use the same format to present 
our results for a given observable in four plots unless 
otherwise mentioned. Our presentation is illustrated in 
Fig. 11 In the two upper panels we present the ratio 
of Eq. ( jll ), as a function of the insertion time (tins), 
shifted by half the sink-source separation, i.e. tins ~ts/2. 
This way, the midpoint time of the ratio coincides for 
all sink-source separations at tins — ts/2 = 0. In what 
follows all times are measured relative to to and thus we 
drop the reference to to. In the third panel we show the 
summed ratio as a function of the sink time, as obtained 


by Eq. (28) and in the bottom panel we compare results 


from the summation method and from the two-state fit 
method with those obtained by the plateau method. 

Let us first discuss the results for the scalar charge 
shown in Fig. H In the two upper panels we show 
the ratio for the isoscalar and isovector scalar charges, 
for the three sink-source separations considered. As ex¬ 
plained in the previous section, when the time separa¬ 
tions Atins ^ 1 and A(G — tins) ^ 1, the ratio becomes 
time-independent. Fitting in the plateau region to a con¬ 
stant value, which we refer to as the plateau value, we 
obtain gs, as in Eq. (27). This is shown by the blue band 
in Eig. Ufor ts = 14a. One observes an increasing trend 
for the plateau value and a clear curvature, especially for 
the isoscalar, indicating dependence on excited states. 
Carrying out a two-state fit yields the dashed lines. As 
in the case of the summation method, the contact points 
tins = to and tins = ts are omitted. The value for gs ob¬ 
tained by the two-state fit is given by the dashed line that 
spans the entire x-range of the figure, with the red band 
indicating the statistical error, while the result of the 
summation method is shown with the solid line and gray 
band indicating the error. As can be seen, within errors 
the two-state fit is consistent with the plateau fit, but 
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FIG. 4: Results for the isovector and isoscalar nucleon scalar 
charge: Upper two panels is the ratio from which gs is ex¬ 
tracted as a function of hna — ts/2 for the isoscalar (upper) 
and the isovector (lower). The blue bands spanning from 
(tins — taji)!a=-4: to 4 are fits to the ratio for ta/a = 14. 
The dashed lines show the result of the two-state fit method. 
The dashed (solid) line spanning the entire x-range show the 
value obtained via the two-state (summation) method, with 
the band indicating the corresponding statistical error. In the 
third panel, the summed ratio is shown for the isovector (filled 
symbols) and isoscalar (open symbols) case. The line shows 
the result of a linear fit, while the bands show the statistical 
error based on the jackknife error of the fitted parameters. In 
the bottom panel, we show the result for gs when using the 
plateau method with ta/a = 10,12 and 14 (squares, circles, 
and rhombuses respectively), as well as when using the sum¬ 
mation method denoted by “sm” (asterisks) and the two-state 
fit “2-st.” (triangles). 


not with the value from the summation method, which, 


FIG. 5: Results for the axial charge. The notation is the same 
as that in Fig. 


however, carries a very large error indicating the need 
to increase statistics in order to have a better assess¬ 
ment of the result of this method. In the third panel of 
Fig.g] we show the summed ratio for the scalar charge 
as a function of the sink time, as obtained by Eq. (281, 
for the isovector and the isoscalar cases. Fitting to a 
linear dependence with respect to ts, one extracts the 
desired matrix element from the slope (Eq. (29)), which 
is the result shown by the gray band in the two upper 
graphs of Fig. The width of the bands is obtained 
by a jack-knife re-sampling of the summed ratio to ob¬ 
tain jack-knife errors for the slope, A4, and intersection 
C of Eq. (29). The values for gs from the summation 
method and those obtained by the plateau method and 
the two-state fit method are shown in the bottom panel 
of the figure. One clearly observes the increasing trend 
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FIG. 6: Results for the tensor charge. The notation is the 
same as that in Fig. 


of the plateau values with increasing t^ja as well as the 
larger values of the summation method shown by the as¬ 
terisks. This study shows that both larger sink-source 
time separations as well as larger statistics are needed 
in order to obtain a meaningful convergence of all meth¬ 
ods. This corroborates our findings of our high statistics 
analysis of the Nj = 2-|-l-l-l TMF ensemble with pion 
mass 373 MeV, referred to as B55.32 ensemble, where 
we showed that ~ 1.5 fm is needed m- Our current 
statistics do not allow to use such a large sink-source 
separation for the physical ensemble. We note that as 
a check of the robustness of the two-state fit, we omit 
more points besides the time slice of the source and the 
sink. In the case of the scalar charge, taking the fit range 
tins G [^0 + 2,ts - 2] and tins G [^o + 3,ts - 3] we ob¬ 
tain: 2.18(34) and 2.22(33) respectively for the isovector 


case and 9.68(26) and 9.75(24) for the isoscalar, which 
are consistent with 2.16(34) and 9.62(27) extracted when 
just omitting the source and the sink. Thus for the scalar 
charge, the fluctuation of the central value when chang¬ 
ing the ht-range is within the statistical error and of the 
order of 2%. 

In Fig. we show results for the axial charge follow¬ 
ing the same notation as that in Fig. For the axial 
charge, one observes a milder dependence on show¬ 
ing that excited states contributions are suppressed for 
this observable. Because of this weaker dependence a 
two-state fit does not yield a meaningful result for these 
values of tg/a, at least within the statistical accuracy of 
1536 measurements. We therefore only show results for 
the plateau and summation methods. The values from 
the plateau method do not vary as a function of tg and 
are in agreement with the value extracted from the sum¬ 
mation method, within the large statistical uncertainties 
of the latter. 

Our results for the tensor charge (gr) are shown in 
Fig.i The dependence of gr on tg is similar to that 
observed in the case of gA and thus a two-state fit fails 
to accurately resolve the excited states. Thus we only 
compare the results extracted using the summation and 
plateau methods in the bottom panel. As can be seen, 
gx exhibits no dependence on the sink-source separation 
within the current statistical uncertainties, evident by 
the same values extracted by fitting to the plateau for 
the three different sink-source separations. The value 
extracted from the summation method is in agreement 
but carries a much larger error and thus does not provide 
a stringent check. 


C. Nucleon momentum fraction, helicity and 
transversity moments 

The momentum fraction is shown in Fig. for the 
connected isoscalar {x)u+d and isovector {x)u-d combi¬ 
nations. Both isoscalar and isovector channels exhibit 
excited state contributions, especially for the isoscalar 
channel. A two-state fit is performed using ting € [to + 
2, tg — 2], which however yields too large errors to include 
in the plots (see Table [Iv| . Using tins G [to -I- 3,ts — 3], 
we find a value of 0.48(19) for the isoscalar, which is con¬ 
sistent with the value given in Table EYl For the isovec¬ 
tor both fit ranges yield consistent results albeit with a 
large error that does not allow to access the sensitivity on 
the fit range. Furthermore, as can be seen in the lowest 
panel of panel of Fig. [^ a decreasing trend is observed 
as the sink-source separation is increased from 10a to 
14a, showing that elimination of excited state effects is 
responsible for reducing the value of this matrix element. 
The summation method yields a value that is even lower 
but with a large statistical uncertainty. Like for the case 
of the scalar charge a larger value of tg/a and increased 
statistics will be needed to reach consistency among the 
various methods with meaningful errors. 
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FIG. 7: Results for the momentum fraction. The notation is FIG. 8: Results for the nucleon helicity moment. The nota- 
the same as that in Fig. tion is the same as that in Fig. 


The helicity moment {x)/^u±Ad and transversity mo¬ 
ment {x)su±Sd are shown in Figs. and For both 
isoscalar observables, a milder dependence on the sink- 
source separation is observed. This is also true for the 
isovector transversity moment. On the other hand, the 
isovector helicity moment shows a decreasing trend sim¬ 
ilar to that observed in the case of the isovector momen¬ 
tum fraction. The value obtained using the summation 
method is consistent in all cases with the plateau value 
when tg/a = 14, albeit with a large statistical uncer¬ 
tainty. 

The nucleon results presented in Figs.|^to|^and Figs.[^ 
to 1^ are summarized in Table [V] where we give the values 
obtained when using the plateau method and the sum¬ 
mation method. Our results with sink-source time sepa¬ 
ration 12a, 14a, and the summation method agree within 


TABLE IV: Results for the nucleon scalar charge and momen¬ 
tum fraction when employing the two-state fit. 

Moment isovector isoscalar 
gs 2.16(34) 9.62(27) 

{x)q 0.19(24) 0.50(20) 


one to two standard deviations. The errors exhibited by 
the summation method, however, are still large, which is 
explained by the fact that this method relies on a two 
parameter fit, contrary to the plateau method which is a 
fit to a constant. Nevertheless, within our current statis¬ 
tics, the summation method can provide an additional 
check of excited state effects. 

We observe that the scalar charge gs and momentum 
fraction {x)u±d exhibit non-negligible excited state effects 
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FIG. 9: Results for the transversity moment. The notation is 
the same as that in Fig. 


when increasing the sink-source separation. For these 
cases we show in Table [TVl the results when employing the 
two-state fit method. For both observables, the two-state 
fit result agrees with the plateau method for tg/a = 14. 

In Table VI we give our results for the isovector quan¬ 
tities as determined from the plateau method using tg 
1.3 fm, as well as for the up- and down-quark contribu¬ 
tions neglecting disconnected diagrams. We note that for 
the up- and down-quark contributions of these quantities 
we carry out the complete analysis with jack-knife resam¬ 
pling starting from three-point correlation functions with 
only an up- or down-quark insertion. Alternatively one 
can form linear combinations of the final isovector and 
isoscalar results of Table |Vj which will give consistent 
up- and down-quark contributions within statistical er¬ 
rors. Except for the scalar and the momentum fraction, 


TABLE V: Results for the nucleon charges and first moments 
computed with the physical ensemble. The first column de¬ 
notes the observable, with u—d indicating the isovector contri¬ 
bution and u -\- d the connected isoscalar. Results extracted 
using the plateau method are given for tg/a = 10,12 and 
tg/a = 14 in the second, third, and fourth columns. The 
value extracted using the summation method is given in the 
last column. The errors are obtained using jack-knife. 


Moment 

10a 

Plateau 

12 a 

14a 

Summation 

Method 


6.46(27) 

7.84(48) 

8.93(86) 

14.0(2.2) 


0.55(18) 

1.18(34) 

2.20(54) 

5.0(1.4) 

9't^ 

0.583(14) 

0.597(23) 

0.611(48) 

0.637(92) 

97^ 

1.158(16) 

1.162(30) 

1.242(57) 

1.28(12) 

9-+'^ 

0.596(21) 

0.598(31) 

0.555(63) 

0.57(12) 

U — d 

9t 

1.062(21) 

1.058(35) 

1.027(62) 

0.99(13) 

{^'ju + d 

0.645(13) 

0.587(18) 

0.540(28) 

0.389(66) 

( 2 ^) u — d 

0.248(09) 

0.218(15) 

0.208(24) 

0.116(54) 

{x)Au+Ad 

0.161(12) 

0.143(17) 

0.175(30) 

0.177(63) 

{x} Au — Ad 

0.286(11) 

0.270(16) 

0.229(33) 

0.159(62) 

Su-\-Sd 

0.195(15) 

0.178(21) 

0.212(36) 

0.204(87) 

{x}su-Sd 

0.311(13) 

0.261(19) 

0.306(29) 

0.245(65) 


TABLE VI: Results for the nucleon axial and tensor charges 
and first moments of parton distributions. In the first column 
we give the observable, in the second the isovector combina¬ 
tion and in the third and fourth the u and d values neglecting 
disconnected contributions. 

isovector up down 

qa 1.242(57) 0.926(47) -0.315(24) 
gr 1.027(62) 0.791(53) -0.236(33) 

(x)q 0.208(24) 0.373(22) 0.166(13) 

{x)Aq 0.229(30) 0.202(26) -0.027(16) 

{x)sq 0.306(29) 0.264(25) -0.045(21) 


the results of Table |VI| are in agreement with the value 
obtained using the summation method. For the scalar 
and the momentum fraction they are consistent with the 
result extracted using the two-state fit albeit with large 
statistical error especially for the momentum fraction. 
Since for the scalar there are large differences still be¬ 
tween the results at different values of tg as well as from 
the value extracted using the summation method we do 
not include a single value in the table. 

Our TMF results for the three isovector charges 
and g^~‘^ and moments {x)u-d, {x)Au-Ad and 
{x)su- 5 d are collected in Figs. 10 and In the Ap¬ 
pendix we give our updated results for our high statis¬ 
tics analysis of the B55.32 ensemble for several sink- 
source separations in Tables |VIII| and |IX| as well as for 
one Nf = 2 -I- 1 + 1 ensemble at a finer lattice spac¬ 
ing [25] in Table These results use the new renor- 

ml For 5 “-^ {x)u-d 


malization functions given in Table 
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FIG. 10: Isovector nucleon scalar charge (upper), axial 
charge (middle), and tensor charge g^~‘^ (lower) using 

the values of Table [Vl] Twisted mass fermion (TMF) results 
are shown for i) Nf = 2, a — 0.089 fm (open green squares), 
a = 0.07 fm (open blue diamonds) and a = 0.056 fm (open 
magenta circles); ii) Nf = 2 + 1 + 1, a = 0.082 fm (filled green 
squares), a = 0.064 fm (filled blue diamonds) and iii) Nf = 
2 TMF clover-improved a = 0.093 fm (physical ensemble), 
ts/a = 12 (filled red triangle), ts/a = 14 (open red triangle), 
summation method (open right triangle). The physical value 
is shown with the black asterisk. For the scalar charge we 
show with the open yellow square the value when ta ~ 1.5 fm. 


and {x)Au-Ad we include results using Nf = 2 at three 
lattice spacings and, for one mass, at two different vol¬ 
umes [351 [2S] • These are given with the new renormal¬ 
ization functions in Table |XI| of the Appendix. These 
results show that cut-off effects are small for lattice spac¬ 
ings smaller than 0.1 fm. Finite volume effects are not 
visible within our statistical accuracy when comparing 
results for two ensembles simulated at a pion mass of 
about 300 MeV and LniTr = 3.3 and Lm^r = 4.6. For the 
physical ensemble we show results for tg = 12a ~ 1.1 fm 



FIG. 11: Isovector nucleon momentum fraction {x)q (upper), 
helicity {x)Aq (middle), and transversity {x)sq (lower). The 
notation is the same as that in Fig. |10| 

and ts = 14a ~ 1.3 fm and from using the summation 
method. We expect {x)u-d to have moderate excited 
states contamination as revealed by our high-statistics 
investigation of {x)u-d for the Nf = 2 + 1 + 1 ensemble 
with m,n. = 373 MeV that showed that its value decreases 
with increasing tg [25] . Our current results at the phys¬ 
ical point with 1536 statistics have much larger errors 
as compared to what was achieved in Ref. [3S] but there 
is a clear decreasing trend as we increase ts- For the 
scalar charge the excited state contributions are large, 
as can be seen both from the results obtained with the 
Nf = 2 + 1 + 1 ensemble with 373 MeV where a high 
statistics analysis is carried out but also for the physi¬ 
cal ensemble where no convergence is achieved with the 
current sink-source separations and statistics. Although 
for {x)u-d and {x) Au-Ad our results from the sum¬ 

mation method are in agreement with the experimental 
values, the errors are still too large and must be reduced 
by a factor of at least two to draw a safe conclusion. How- 
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TABLE VII: Results for the renormalized in MS at 

2 GeV for the four TMF clover-improved ensembles consid¬ 
ered in this work. In the first column we give the bare twisted 
quark mass, and in the second column the spacial extent in 
lattice units. In the fourth column we provide the value of 
{x)1,-d of the pion with its statistical and systematic uncer¬ 
tainty computed as explained in the text. The last column 
gives the number of measurements for each ensemble. 


a/L 

L/a 

{x}Z-d in MS at 2 GeV Ameas 

0.006 

24 

0.249(7)(133) 

210 

0.006 

32 

0.259(4)(t() 

240 

0.003 

24 

0.219(13)(/®) 

276 

0.0009 

48 

0.214(15)(t(,") 

309 



ever, we stress that our value for from the plateau 
method using ts ~ 1.3 fm agrees with the experimental 
value. To our knowledge, this is the first computation for 
which the value of the axial charge is reproduced from the 
plateau method, without any chiral extrapolation. 


D. Pion momentum fraction 


In this section we present results on the isovector pion 
momentum fraction. Three Nf = 2 TMF ensembles with 
the clover term are analyzed with heavier than physical 
pion masses, two of which with spatial lattice size 2.23 fm 
and one with spatial lattice size 2.98 fm. A fourth en¬ 
semble that includes the clover term is simulated using 
the physical value of the pion mass and spatial lattice 
extent of 4.46 fm. This is the ensemble used for the nu¬ 
cleon observables and the ensemble details can be found 
in Table |T] The number of measurements, which are well 
separated in the number of HMC traiectories, is given in 
Table IvTtI 


In Fig. 12 we show the ratio Eq. (22) as a function of 


tins/a for the physical ensemble. The black horizontal 
line represents the value quoted in Table VII The sta¬ 
tistical accuracy of the pion correlation functions allows 
for a more careful assessment of systematic uncertainties 
in the plateau fit as compared to the case of the nucleon. 
Namely, we obtain the plateau value by performing con¬ 
stant fits to the data with all possible fit ranges with 
degrees of freedom larger than 5. For each of these fits a 
weight factor 


l-2|p-0.5|y /I-2|p„„ -0.5P " 


W 




is computed, where p (pm„) is the p-value of the fit and 
W (Wm^) the statistical error of the fit parameter (of 
rriTr) determined using 1500 bootstrap samples. The pion 
mass itself is also determined for a large number of fit- 
ranges. The final result is determined as the weighted 
median over all combinations of fit-ranges. The 68.54% 


FIG. 12: The ratio for {x)'^d using the ap = 0.0009 ensemble. 
We show the weighted median over the different ht-ranges as 
a solid black line, the statistical error as the red band and the 
systematic errors as the gray band. 



suits obtained for lattices of two sizes: L = 24a (red squares) 
and L = 32a (blue circles). 


confidence interval of the weighted distribution is quoted 
as the systematic uncertainty. 

For the twisted mass value a/i = 0.006 we have two 
spatial lattice sizes available, namely L/a = 24 and 
L/a = 32. Within errors, the result for agrees 

between these two ensembles. We show the re¬ 

sults for L/a = 24 and L/a = 32 in Fig. 13 as a function 

of tins/ 

We compare our Nf = 2 clover-improved results with 
results obtained for Nf = 2 twisted mass ensembles with¬ 
out the clover term published in Ref. (SS]. These ensem¬ 
bles were simulated using /3 = 3.90 with four values of the 
bare quark mass: ap = 0.004, 0.0064, 0.0085 and 0.0100. 
The lattice spacing for these ensembles is 0.089(1) fm, 
which is similar to the clover-improved ensembles. For 
a/i = 0.004 we have again two spatial lattice extents 
L/a = 24 and L/a = 32 available. As for the clover- 
improved ensembles the results for agree between 

the two volumes. 
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2 GeV using three loop perturbation theory. The phe¬ 
nomenological value is compatible with the value of 
0.214(15)(tg^) computed for the physical ensemble. 

The results presented here can be compared to Ref. [5] , 
where Nf = 2 non-perturbatively clover improved Wilson 
fermions have been used, including two ensembles with 
pion mass values around 150 MeV. Two values of the lat¬ 
tice spacing are investigated, a = 0.06 and a = 0.07 fm, 
respectivley. In that reference a bending of the momen¬ 
tum fraction values towards small pion mass values is 
observed, while the agreement at > 300 MeV to the 
results presented here is reasonable. 


E. Comparison of nucleon observables with 
different fermion actions 


FIG. 14: The renormalised momentum fraction of the pion 
{x)Z-d a function of the squared pion mass at renormal¬ 
ization scale p = 2 GeV in the MS scheme. The results of 
this work (red symbols) are shown together with previous re¬ 
sults obtained using Nf = 2 ensembles (green symbols) |36| 
compared to the phenomenological value (black star) from 
Ref. |64] . Results at two values of the pion mass but different 
lattice volumes are shown by the open squares (24^ x 48) and 
diamonds (32^ x 64). 


From the comparison of the different available spatial 
lattice sizes we conclude that within the current statisti¬ 
cal uncertainties we cannot detect significant finite vol¬ 
ume effects for rriTrL < 3.2 realized for the /3 = 3.90, afi = 
0.004, L/a = 24 ensemble. The physical ensemble has 
slightly smaller niTrL = 2.97, and the clover ensemble 
with afx = 0.003 has = 2.77. Therefore, we can¬ 
not completely exclude finite size effects for these two 
ensembles. All other ensembles have rriTrL > 3.2. Note 
that we are currently generating a physical ensemble with 
Lja = 64, which will allow us to check for finite size ef¬ 
fects. We expect pion observables to be more sensitive 
to finite size effects than nucleon observables; having a 
larger volume will enable us to confirm this expectation. 

The results for the renormalized isovector mome utun 
fraction of the pion are summarized in Table VII 

The renormalization factors used are given in Table m 
at 2 GeV in the MS scheme. The results are also dis¬ 
played in Fig. as a function of the squared pion mass. 
In Fig, one observes that there is agreement between 
the clover-improved and non-clover improved ensembles 
within errors. Note that systematic uncertainties are not 
displayed. 

In Fig. we compare with the latest phenomenologi¬ 
cal value for which can be found in Ref. [64] and 

reads 

= 0.256(13). 

Note that the result given in Ref. [M] is at /r = 5.2 GeV 
renormalization scale and we have translated it to /x = 


In this section we compare our results on nucleon ob¬ 
servables with other recent results obtained using simu¬ 
lations with similar parameters. 



FIG. 15: Results for the nucleon axial charge for different 
fermion actions. Twisted mass fermion results are shown with 
open green squares for Nf = 2 ensembles |22| . with filled 
blue squares for Nf = 24-1-1-1 [26] and with the open red 
triangle for the physical ensemble using the plateau value at 
ts/a — 14 (see Table VI I. Results are also shown using Nf = 2 
clover fermions (filled purple diamonds) |65| : Nf = 2 4-14- 
1 staggered sea and clover valence quarks (Hlled light blue 
inverted triangles) (SHI ; Nf = 24-1 with DWF on a staggered 
sea (Hlled yellow circles) [63; and Nf =24-1 clover (black 
x-symbols) [68]. 


A number of lattice QCD collaborations are investigat¬ 
ing gA since, as emphasized already, this is considered a 
benchmark quantity for lattice QCD. In Fig. [^we show 
results for Nf = 2 [32] and V/ = 24-14-1 [23 twisted 
mass fermions obtained in previous analyses using sim¬ 
ulations with pion masses in the range of 450 MeV to 
210 MeV for various volumes always satisfying the con¬ 
dition > 3. For Nf = 2 ensembles, three values of 
the lattice spacing were analyzed, namely a = 0.089 fm, 
0.07 fm and 0.056 fm and, at one pion mass of about 
300 MeV, for two different volumes. As already pointed 
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FIG. 16: Isovector nucleon scalar charge (upper) and 
tensor charge (lower) versus m^. Twisted mass fermion 
(TMF) results are shown for two ensembles ofA^f = 2-|-l-|-l 
fermions (filled green square for ta ~ 1.15 fm and open green 
square for ta ~ 1.5 fm, and filled blue diamond) and for the 
physical ensemble (filled red triangle for ta ~ 1.1 fm and open 
red triangle for ta ~ 1.3 fm). Results are also shown using: 
clover fermions on Nf = 2 + 1 -|- 1 staggered sea from Ref. |66l 
169) for gg~‘^ and from Ref. [70], for g^~'^ (filled light blue 
downwards triangles); Nf =2-1-1 clover (black filled triangles 
and crosses), Nf = 2-1-1 domain wall fermions (open light blue 
circles) and hybrid (blue crosses) |7T] ; Nf = 2 Clover fermions 
for three values of the lattice spacing (filled magenta, yellow 
and light blue circles crosses) |65) . All results were extracted 
using the plateau method except those from Ref. [Ml [70], 
which used a two-state fit. 


FIG. 17: Isovector nucleon momentum fraction {x)n-d and 
helicity {x}Au-Ad- Twisted mass fermion results are shown 
for Nf = 2 ensembles (open green squares), for two Nf = 
2 -I- 1 -I- 1 ensembles (blue filled square) and for the phys¬ 
ical ensemble with a clover term (open red triangle) taken 
from Table |VI| Also shown are results from RBC-UKQCD 
using Nf = 2-1-1 DWF (magenta right pointing triangle) | 72 |. 
from LHPC using DWF on Nf =2-1-1 staggered sea (blue 
crosses) |67| and QCDSF/UKQCD using Nf = 2 clover 
fermions (filled magenta diamond) [73]. For (x)u-d we also 
show results from LHPC using Nf =2-1-1 clover with 2-HEX 
smearing (filled black triangles) [68] and Nf = 2 clover (open 
black circle) |74| . All values are extracted using the plateau 
method and ts ~ (1 — 1.2) fm, except our result at the phys¬ 
ical point for which G ~ 1.3 fm was used. The experimental 
value for {x)u-d is taken from Ref. [75] and for {x)Au-Ad from 
Ref. [76]. 


out, the consistency among these results indicates small 
cut-off and finite volume effects. The Nf = 2 val¬ 
ues are consistent with the values extracted using two 
Nf = 2-I-1-I-1 ensembles with lattice spacing a = 0.082 fm 
and a = 0.064 fm, showing that there are no visible 
strange and charm sea quark effects on these quantities 
at least to the accuracy we now have. This allows a com¬ 
parison with results using different fermion discretiza¬ 
tion schemes even before the continuum extrapolation is 
performed. In Fig. we include results obtained us¬ 
ing clover improved fermions from two collaborations: 
In Ref. [S^ results were obtained using Nf = 2 clover 


fermions with smallest pion mass of 150 MeV and three 
lattice spacings a = 0.08 fm, 0.07 fm and 0.06 fm as well 
as several volumes. These results are in agreement with 
ours. The LHPC analyzed Nf = 2 + \ tree-level clover- 
improved Wilson fermions with 2-HEX stout smeared 
gauge links provided by the BMW collaboration using 
smallest pion mass of 149 MeV at one lattice spacing 
a = 0.116 fm [68]. These tend in general to have lower 
values. This is particularly severe close to the physical 
pion mass. LHPC also computed the axial charge in a 
mixed action approach that uses DWF on staggered sea 
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quarks by LHPC m with a = 0.124 fm where high ac¬ 
curacy results were produced for heavier pion masses. 
These results are in good agreement with ours. We note 
that both TMF and clover-improved results are extracted 
using the plateau method with sink-source time sepa¬ 
ration of about 1 fm to 1.2 fm with the exception of 
our result for the physical ensemble where we used time 
separations of up to 1.3 fm. In addition, results at two 
pion masses using a hybrid action with clover valence on 
Nf = 2-I-1-I-1 staggered fermions used a two-state fit [66] 
and are included here for comparison. They tend to be 
higher than other results although they are compatible 
with our Nf = 2-I-1-I-1 value at pion mass of about 
210 MeV, which albeit carries a large error. The general 
conclusion is that the lattice results for pion mass higher 
than about 300 MeV, which are more accurate as com¬ 
pared to those at smaller pion masses, are in agreement. 
This is an indication that lattice systematics are under 
control in this pion mass range. Results for pion masses 
smaller than about 300 MeV, in general, have larger sta¬ 
tistical errors and agreement among them is harder to as¬ 
sess. They clearly indicate the need for more precise val¬ 
ues and a reliable assessment of systematic uncertainties. 
This is particularly relevant for the results close to the 
physical point where we observe a disagreement between 
clover results from LHPC at pion mass of 149 MeV and 
from Ref. [65| at similar pion mass. This discrepancy was 
claimed to be due to excited states, which were shown to 
be suppressed with improved smearing in Ref. [65j . This 
needs to be further investigated with a dedicated preci¬ 
sion calculation with a complete assessment of system¬ 
atic uncertainties. Other results using clover-improved 
fermions not shown here are those by the CLS collabora¬ 
tion |40j . which extracted their values from the summa¬ 
tion method. A complete set of the results on gA can be 
found in Ref. |4]. The result of this work is shown with 
the open triangle obtained for tg = 1.3 fm. This value is 
in agreement with the experimental value with, however 
admittedly a rather large error. 


The calculation of the scalar and tensor charges has 
received more attention recently due to their relevance 
for searches of new scalar and tensor couplings beyond 
the familiar weak interactions of the Standard Model in 
the decay of ultra-cold neutrons. We compare our TMF 
results in Fig. with those obtained by three groups 
whose results on the nucleon axial charge were also in¬ 
cluded in Fig. 15 The first set of results are from Ref. [HS] 
using Nf = 2 clover fermions at three lattice spacings. 
The second set is from the LHPC group which in Ref. m 
used V/ = 2 -I- 1 clover with 2-HEX stout smeared gauge 
links at lattice spacings a = 0.116 fm and a = 0.09 fm, 
Nf = 2-1-1 DWF with a = 0.084 fm and a hybrid ac¬ 
tion of DWF on staggered sea with a = 0.124 fm. Both 
these groups used the plateau method and sink-source 
time separation within 1 to 1.2 fm. The third set of 
results are from Ref. [66] at one lattice spacing for the 
scalar and from Ref. m at three lattice spacings for 
the tensor obtained using a hybrid action of DWF on 


7Vy = 2-|-l-|-l staggered fermions and employing a two- 
state fit. For the case of the scalar charge shown in the 
upper panel of Fig. [^ we observe overall a good agree¬ 
ment among all lattice results obtained with similar sink- 
source separation. However, our high-statistics analysis 
using Nf = 2-I-1-I-1 TMF at pion mass 373 MeV revealed 
excited states contamination, which only become negli¬ 
gible when tg ~ 1.5 fm, increasing the value of gs- The 


value extracted when tg = 1.48 fm is shown in Fig. 16 


by the open green square. A similar analysis for the 
physical ensemble also reveals large contributions from 
excited states for gg~‘^. Comparing results obtained for 
tg ~ 1.1 fm and 1.3 fm we confirm an increasing value 
as we increase tg. However, the statistical error is large 
despite the fact that we have 1536 measurements as com¬ 
pared to 1200 used for the ensemble at ttItt = 373 MeV. 
This demonstrates that, obtaining the same accuracy at 
the physical point for tg ~ 1.5 fm, which may be needed 
to suppress excited states, requires more than an order 
of magnitude increase in statistics. 

Results on the isovector tensor charge are compared 
in the lower panel of Fig 16 Our TMF results shown 
in Fig, show that excited state contributions are less 
severe for g^~‘^ and that the values at tg/a = 12 and 
tg/a = 14 are consistent. Indeed our value at the phys¬ 
ical point obtained using tg ~ 1.3 fm is in very good 
agreement with other lattice results providing a predic¬ 
tion for this important quantity directly at the physical 
point. 


Recent lattice QCD results have also been obtained for 
the isovector momentum fraction and helicity. A compar¬ 
ison of our results for {x)u-d and {x)Au-Ad with other 
collaborations is shown in Fig. [IT] We only show results 
extracted using the plateau method. Most of the analyses 
employed a sink-source separation of 1 to 1.2 fm includ¬ 
ing our TMF Nf = 2 ensembles. As shown in Ref. [25] 
where {x)u-d was computed using our Nf = 2 -|- 1 -f 1 
ensemble at pion mass of 373 MeV and high-statistics, 
excited states may not be negligible for this observable. 
Indeed, this is confirmed by our current analysis for the 
physical ensemble where there is a decreasing trend as 
tg increases. Our value at the physical point is in agree¬ 
ment with the other lattice values extracted close to the 
physical point. These are from Ref. [74], which it is an 
update of Ref. [77] for 160 MeV and from LHPC 

at TOtt 150 MeV using Nf = 2-1-1 clover fermions 
with 2-HEX smeared gauge action [^. Our value at 
tg ^ 1.3 fm is still larger than the experimental value. 
We are currently performing a high statistics analysis for 
our physical ensemble using larger values of tg to inves¬ 
tigate contamination due to excited states, which tend 
to decrease the value of {x)u-d- For larger pion masses 
we show results using Nf =2-1-1 DWF from the RBC- 
UKQCD collaborations [75], from LHPC [ST] using DWF 
on Nf = 2 + 1 staggered sea and from the QCDSF collab¬ 
oration using Nf = 2 clover fermions [73]. Results from 
LHPC used perturbative renormalization which could ex¬ 
plain the fact that these are in general lower than other 
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lattice results. For the case of {x)Au-Ad the situation is 
similar and our value using the physical ensemble is still 
larger than its experimental value. As for the momentum 
fraction there is a decreasing trend as ts increases. In fact 
the summation method yields a value that is consistent 
with the experimental value as can be seen in Fig. [m 
However, the error is too large and our goal in a future 
analysis is to reduce it by a factor of two so as to confirm 
agreement with the experimental value. Resolving these 
discrepancies will give more confidence on our prediction 
for the transversity moment. 


VI. CONCLUSIONS 


In this work we present results on the pion momentum 
fraction and key nucleon observables using lattice QCD 
simulations at the physical value of the up and down 
quarks. Our analysis of the isovector pion momentum 
fraction uses Nf = 2 ensembles with the clover term sim¬ 
ulated at three different values of the light quark mass. 
We find a value of {x)^^ = 0.214(15)(t9^) in the MS at 
2 GeV at the physical point. 

For the nucleon system, we compute the three local 
and three one-derivative isovector and isoscalar matrix 
elements at zero momentum transfer. In our calculation 
we analyze three sink-source time separations, which al¬ 
lows us to investigate excited state effects by observing 
the dependence of the extracted nucleon matrix elements 
on this separation. For all observables we compare the 
plateau method with the summation method. In some 
cases the sensitivity on the sink time is good enough so 
that a two-state fit can also be applied as a third method 
to detect excited state contaminations. Employing these 
different methods to ensure that contamination from ex¬ 
cited states is suppressed is crucial in obtaining reliable 
results. However, for this study to be meaningful one 
needs large statistics in particular for large sink-source 
time separations and for the summation method. For 
the pion momentum fraction where statistical errors are 
smaller one extracts the relevant matrix element using 
the largest possible time separation ensuring ground state 
dominance. Our results for the nucleon axial charge and 
isovector pion momentum fraction are in agreement with 
their experimental values, that constitutes a very impor¬ 
tant conclusion of this study. Since the tensor charge is 
found to behave similarly to the axial charge as far as 
ground state dominance is concerned we can predict its 
value at the physical point to be g^~‘^ = 1.027(62) in 
the MS scheme at 2 GeV. Assuming that disconnected 
contributions remain as small at the physical point as 
were found at a pion mass of 373 MeV where they were 
shown to be negligible EH EH], we can give a direct 
prediction for the individual up- and down-quark tensor 
charges. We find = 0.791(53) and g^ = —0.236(33) 
(see Table VI). 

Thus, this first lattice study of nucleon and pion struc¬ 
ture at the physical values of the light quark masses is 


very promising for future precision calculations of these 
key quantities directly at the physical point. Ongoing 
plans include an analysis with increased statistics for gen¬ 
eral momentum transfer and new 7V^ = 2-|-1-I-1 simula¬ 
tions with their mass fixed to physical values, combined 
with larger volumes and improved algorithms for noise 
reduction such as multiple right-hand-side solvers. Af¬ 
ter reproduction of benchmark quantities such as gA for 
the nucleon and the pion, lattice QGD is in a position 
to turn to quantities more difficult to obtain experimen¬ 
tally such as the scalar and tensor charges gs and gr- 
Such charges are of high interest in phenomenology and 
experiments since these enter in couplings of protons to 
super-symmetric candidate particles. Their precise de¬ 
termination can therefore be used to exclude regions in 
dark matter searches and influence future experimental 
set-ups for new physics searches. 
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Appendix 

In this Appendix we give the ETMC results for the 
Nf = 2 and the Nf = 2 + 1 + 1 ensembles published in 
Refs. ESI respectively. These results are updated us¬ 
ing the new renormalization functions given in Table |lIT| 
In Tables lyHIl and HXl we collect the results for the B55.32 
ensemble for which a high statistics analysis is carried out 
for several sink-source separations. In Figs. [TS] and [T^ we 
show the ratios from which the isovector and isoscalar 
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connected charges and first moments of PDFs are ex¬ 
tracted. The ratios are computed for a total of five sink- 
source time separations spanning a time range of about 
0.8 fm to 1.5 fm enabling us to apply the summation 
method as check for . This high statistics analysis allows 
us to perform a two-state fit for all quantities except the 
axial charge where the excited state contamination is the 


mildest. As our final values we take the plateau value 
that is in agreement with the value extracted from the 
summation method and two-state fit when possible. Fi¬ 
nally in Tables and |X1| we give the results for all our 
other ensembles where only one sink-source separation 
was employed. 


TABLE VllI: Updated results for the Nf = 2 + 1 -|- 1 B55.32 ensemble for the nucleon charges and first moments of parton 
distribution functions. For the isoscalar combination we give only the connected contribution. 


ts/a : 10 

Plateau method 

12 14 

16 

18 

summation 

method 

two-state 

fit 

isovector 1.08(3) 
isoscalar 5.07(4) 
stat. 2429 

1.12(3) 

5.45(4) 

4396 

1.16(4) 

5.74(6) 

4396 

gs 

1.08(10) 

5.95(13) 

2018 

1.46(20) 

6.33(28) 

1200 

1.19(10) 

6.46(15) 

1.23(10) 

6.81(23) 

isovector 1.143(4) 1.152(5) 
isoscalar 0.596(3) 0.591(4) 
stat. 2429 4396 

1.155(8) 

0.589(7) 

4396 

9A 

1.174(20) 

0.605(16) 

2018 


1.184(19) 

0.583(16) 


qt 

isovector 1.119(6) 1.087(7) 1.058(11) 1.080(30) 
isoscalar 0.680(5) 0.666(6) 0.660(9) 0.663(21) 
stat. 2278 4040 4040 1762 


1.023(27) 

0.624(22) 

1.053(21) 

0.646(9) 


TABLE IX: Updated results for the Nf = 2-|-l-|-l B55.32 ensemble for the nucleon first moments of parton distributions. For 
the isoscalar combination we give only the connected contribution. 


Plateau method 


summation 

two-state 

ts/a : 10 12 14 16 

18 

method 

fit 

{X}q 

isovector 0.290(4) 0.270(3) 0.252(4) 0.233(9) 
isoscalar 0.720(8) 0.677(5) 0.639(6) 0.607(11) 
stat. 2429 4396 4396 2018 

0.252(19) 

0.616(21) 

1200 

0.223(9) 

0.554(15) 

0.214(13) 

0.558(19) 

{x)Aq 

isovector 0.328(3) 0.312(3) 0.297(3) 0.298(8) 
isoscalar 0.207(3) 0.198(2) 0.189(3) 0.193(8) 
stat. 2429 4396 4396 2018 


0.270(8) 

0.172(7) 

0.286(9) 

0.184(7) 

{x}sq 

isovector 0.372(5) 0.349(4) 0.322(5) 0.316(12) 
isoscalar 0.254(4) 0.239(4) 0.219(6) 0.215(15) 
stat. 2278 4041 4041 1763 


0.283(14) 

0.178(14) 

0.284(17) 

0.183(21) 














s^kns/ V.Ls,Lins. 


22 



CO 




(tins-ts/2)/a 


(tins-ts/2)/a 


FIG. 18: The ratios from which the isovector charges (left) and the first moments of PDFs (right) are extracted as a function 
of the insertion-source time separation for the B55.32 ensemble. The statistics used are given in Tables [VIIIl and |IX| 
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FIG. 19: The ratios from which the isoscalar charges (left) and the first moments (right) are extracted as a function of the 
insertion-source time separation for the B55.32 ensemble. Only connected contributions are included. The statistics used are 
given in Tables [Vml and[IXl 
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TABLE X: Updated results for the Nf = 2 + 1 + 1 TMF ensemble with nin = 213 MeV and a — 0.064 fm. In the last column 
we give the number of measurements for all observables. 


/3 {L^ X T) (GeV) 




(^) u — d 

(x) Alt —Ad 

(^) Su—Sd 

ts/a Statistics 

A/=2 + l + l 

2.10 (48^x96) 0.2134(6) 

1.185(61) 

1.024(368) 

1.110(61) 

0.241(19) 

0.276(21) 

0.275(26) 

18 900 


TABLE XI: Updated results for the Nf = 2 TMF ensembles. 


13 (L^ X T) 

(GeV) 


(^) U — d 

Alt— Ad 

ts/a 

Statistics 

Nf=2 

3.90 (24®x48) 

0.3032(16) 

1.129(34) 

0.270(17) 

0.291(10) 

12 

943 


0.3770(9) 

1.158(28) 

0.273(11) 

0.303(9) 


553 


0.4319(12) 

1.152(25) 

0.272(10) 

0.311(8) 


365 


0.4675(12) 

1.181(19) 

0.272(9) 

0.314(6) 


477 

3.90 (32®x64) 

0.2600(9) 

1.174(48) 

0.279(13) 

0.282(14) 


667 


0.2978(6) 

1.121(32) 

0.258(16) 

0.294(10) 


351 

4.05 (32®x64) 

0.2925(18) 

1.211(67) 

0.241(25) 

0.317(24) 

16 

447 


0.4035(18) 

1.191(32) 

0.255(12) 

0.319(10) 


326 


0.4653(15) 

1.189(25) 

0.269(9) 

0.310(7) 


419 

4.20 (32®x64) 

0.4698(18) 

1.138(25) 

0.250(13) 

0.294(9) 

18 

357 

4.20 (48®x96) 

0.2622(11) 

1.142(44) 

0.274(20) 

0.290(15) 
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